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The classical problem of foam film rupture dynamics has been investigated when surfaces 
exhibit very high rigidity due to the presence of specific surfactants. Two new features 
are reported. First a strong deviation to the well-known Taylor-Culick law is observed. 
Then, crack-like patterns can be visualized in the film; these patterns are shown to appear 
at a well defined deformation. The key role of surface active material on these features 
is quantitatively investigated, pointing the importance of surface elasticity to describe 
these fast dynamical processes, and thus providing an alternative tool to characterize 
surface elasticity in conditions extremely far from equilibrium. The origin of the cracks 
and their consequences on film rupturing dynamics are also discussed. 


1. Introduction 

Despite its apparent useless character and simplicity, the dynamics of bursting of soap 
bubbles have fascinated scientists for more than a century. Lucien Bull (1904) made the 
first images of soap bubble bursts. The first theoretical analysis dates back to Dupre and 
then to Taylor (1959) and Culick (1960) where they considered the presence of a rim at 
the edge of a hole created in the liquid film, collecting the liquid during its movement. 
The constant hole opening velocity Vc results from a balance between the rim inertia and 
surface tension in the film, and is given by Vc = with 7eq the equilibrium 

surface tension, p the liquid density and ho the film thickness. These results are in good 
agreement with stationary experiments performed on liquid sheet (Taylor 1959) and has 
been extensively investigated by McEntee & Mysels (1969) in the case of soap films 
thicker than 50 nm. More recently, satellite formation during edge retraction (Lhuissier 
& Villermaux 20096) and bubble entrapment (Bird et al 2010) have been investigated 
as these behaviors are crucial in many applications. Destabilization of liquid sheets or 
bubbles indeed arise in many practical situations ranging from the building material 
industry, when glass sheets are molded, to foam engineering, food processing, biological 
membrane and environmental science (Bird et al 2010). In these applications, liquids 
can be viscous or contain surface active materials. In the latter, surface tension becomes 
a dynamical quantity, which depends on the local surface concentration of surfactants, 
and thus on the elongation of the surface; this is characterized by the surface elasticity 
defined as the derivative of surface tension with respect to relative changes in surface area. 
The effect of surface elasticity has been observed through the development of an aureole 
surrounding the opening hole and expanding with time (Florence & Frens 1972; Liang 
et al 1996; Lhuissier & Villermaux 2009a). However, except in the case of very viscous 
liquid, the opening dynamics always obey Taylor-Culick law, although some deviations 
have been reported by Mysels (McEntee & Mysels 1969; Elorence & Mysels 1974), but 
hardly commented. In this work, we investigate the dynamics of bursting of circular foam 
films generated from surfactant solutions inducing large surface elasticities and we report 
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Figure 1. (a) Image sequence of a foam film rupture (/io=10 /xm, solution E - table 1). The 
timelapse between images is 8 ms. (b) Radius r of the hole vs time t (/io=10 /xm, solution E). 
The red (light) line represents the prediction of Taylor-Culick, while the black line shows the 
initial opening at constant velocity uq. Inset: Picture of a ruptured foam film (solution C) which 
highlights aureole formation, (c) Picture of a foam film (solution E) 37 ms after its breaking. 
The arrows highlight crack-like patterns, which appear during the hole opening. 


for the first time systematic deviations to Taylor-Culick law. A careful analysis allows 
us to estimate surface elasticity at both large compression and compression rate in good 
agreement with reported data in the literature. Moreover, unexpected effects of frame 
size are observed through the appearance of new patterns, reminiscent of fractures or 
wrinkles in the film. 


2. Experimental set-up 

The experimental set-up consists in a circular metallic frame of radius R = 1.5 — 11 cm 
pulled out from a surfactant solution at different velocities to generate films with various 
thicknesses. The film absolute thickness is determined through an absorption technique 
measurement (Lastakowski et al. 2014; Petit et al 2015) and we denote ho the initial 
average thickness of the film. Film rupture is initiated by approaching a heated needle 
and is recorded via a high-speed camera (10000 Hz, Photron SA-4). An image sequence 
is reported in figure la, where we measure the radius r of the expanding hole versus 
time, as shown in figure lb. Surfactant solutions are produced in a 10%-90% glycerol- 
water mixture in which a dye (Brilliant Black BN 60%, Sigma, 5g/L) is added. They 
contain 3.3 g.L“^ of sodium lauryl-dioxyethylene sulfate (SEES, Stepan), 1.7 g.L“^ of 
cocoamidopropyl betaine (CAPB, Goldschmidt) and myristic acid (MAc, Fluka) in the 
concentrations C described in table 1. The surface elasticities of similar solutions are 
well characterized in the literature (Mitrinova et al. 2013a) and span over two orders 
of magnitude when the concentration C of MAc is varied as reported in table 1. Such 
elastic moduli are attributed to the surface properties of the adsorbed layer of MAc, 
whose surface concentration is expected to increase with C up to the saturation of the 
surface (Golemanov et al. 2008). At the same time, micelles of the two co-surfactants 
(SEES and GAPB) help to solubilize the poorly soluble fatty acid. 
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Solution 

A 

B 

C 

D 

E 

F 

C [mM] 

- 

0.055 

0.11 

0.22 

0.88 

2.2 

7eq [mN/m] 

29 

29 

27 

26 

23 

22 

Eod [mN/m] 

4 

50 

90 

200 

400 

400 

Eo(uo) [mN/m] 

60 

200 

2000 

5000 

2.10^ 

4.10^ 

Fo (cracks) [mN/m] 

- 

- 

90 

200 

300 

300 


Table 1. Properties of the surfactant solutions used in the experiments: MAc concentration C, 
equilibrium surface tension 7 eq and surface elasticities. Data from Mitrinova et al (2013a) for 
similar solutions (without glycerol and dye) are reported for 7 eq and the elastic modulus ^od 
measured with the oscillating drop method for small deformation (0.2 —4 %) at frequency 0.2 Hz. 
Eo(uo) corresponds to the elasticity deduced from the initial hole velocity using equation (4.1). 
^0 (cracks) corresponds to the elasticity deduced from cracking radius using equation (4.2). 


3. Results 

Some remarkable features can be underlined. At first, the opening velocity is constant 
as expected but smaller than predicted by Taylor-Culick law (figure lb). Moreover, 
an aureole already described in the past (Florence & Frens 1972; Liang et al. 1996; 
Lhuissier & Villermaux 2009a) is observed through spatial variations of transmitted 
light, especially for the less rigid interfaces (inset of figure lb). Then, some dark patterns 
are observed (see arrows in figure Ic), which we denote cracks in the following. This 
apparition coincides with a decrease of the opening velocity (figure lb), the presence of 
these cracks modifying the bursting dynamics. 

The initial opening velocity uq is represented in figure 2a as a function of the initial 
film thickness ho for various solutions. Without MAc (solution A of table 1), the 
velocity follows Taylor-Culick law (o), which is consistent with interfaces of low elasticity. 
However, in the presence of MAc, the initial velocity is lower than in the previous case. For 
each MAc concentration, the initial velocity varies with l/\/ho oc W- For each solution 
and different thicknesses, we thus extract the initial opening velocity normalized by Culick 
velocity. This quantity decreases when the MAc concentration increases (figure 2b), that 
is, for larger surface elastic moduli (Mitrinova et al. 2013a). 

During the film opening, orthoradial cracks (perpendicular to the direction of opening) 
appear in the film (figures la and Ic), at a well defined radius of the hole r^. Some specific 
irregular fold-like patterns and filaments have been previously reported by McEntee 
& Mysels (1969), although not directly comparable to our observations. For a given 
solution, figure 3a shows that the ratio rp/R is independent of the frame radius (for 
R = l.b — 11 cm) and almost independent of the film thickness (for /iq = 2 — 20 jam). 
The cracks thus appear for a well-defined critical compression of the interface. Figure 3b 
shows that this critical compression decreases with MAc concentration and the surface 
modulus. 


4. Discussion 

These two observations concerning the initial opening velocity and the onset com¬ 
pression for cracks can be rationalized following the framework initially proposed by 
Frankel & Mysels (1969) for the theoretical description of aureoles. They considered 
that surfactants are insoluble, which is reasonable at the timescale considered here: the 
duration of the opening R/uq, typically 30 ms, is smaller than surfactant desorption time 
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Figure 2. (a) Initial opening velocity uq of the hole as a function of the film thickness ho for 
R = 3 cm. The MAc concentration C decreases from dark to light points: solutions A (o), B 
(□), C (>), D (o), E(A) and F (★). (b) initial opening velocity normalized by Culick velocity 
uojVc as a function of C (error bars: 95 % confidence intervals). The magenta solid line shows 
the value measured for C = 0, with error bars shown by the dotted lines. The values uolVc 
are extracted by performing least square percentage fit for each solution, with weights taking 
into account the 1 fim error on thickness measurements. In both figures, the black dashed lines 
represent Taylor-Culick law. 



b) 



Figure 3. (a) Critical compression of the interface for crack formation rpjR as a function of 
the initial thickness ho for solution E and different frame radii = 1.5 cm (o), = 3 cm (□) et 

R = 11 cm (A), (b) rpjR averaged for thicknesses ho = 2 — 35 /xm and R = 3 cm, as a function 
of MAc concentration C. 


r. Indeed, although these processes are likely to be dominated by surfactant exchange 
with micelles in our systems (Golemanov et al 2008), a lower bound for r is provided by 
the diffusion time across the film thickness h^/D ^ 40 ms — 2s (for ho = 2 — 40 /im and 
D = 10“^^ m^/s). Adsorption times longer than 30 ms for myristic acid in these systems 
have also been reported (Mitrinova et al 20135). A compressive shock thus propagates 
at the surface of the film. The liquid is collected in an extended rim — an aureole — 
visible in figure lb (inset) and whose shape depends on surface tension, film thickness 
and surface elasticity. 

Viscous effects have also been neglected. Indeed, as no shear takes place within the 
film thickness, the characteristic Reynolds number and surfaee Reynolds numbers read 
Re = uoRju ^ 1 and Res = puoRho/hz^ respectively, with u the kinematic bulk viscosity 
and K the intrinsic surface viscosity. Surface viscous dissipation can a priori not be 
neglected if values of k measured at 0.2 Hz are considered (Costa 2012; Golemanov 
et al. 2008). However, surface viscosity is expected to collapse at large frequencies, as 
shown in experiments and modeling (Lucassen & Van Den Tempel 1972). Eventually, the 
observation of a constant initial velocity varying with 1/v^ oc Vc (figures lb and 2) is 
a key indication that inertia (and not viscous effects) is dominant in this problem. 
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Figure 4. Solid black line: numerical prediction for the normalized hole velocity 
uqIVc = /(F^o/7eq) for radial bursting (equation (4.1) and A). The dashed line corresponds 
to uojVc = ^/yeq/F^o oxpoctod for unidimensional bursting (Frankel & Mysels 1969). Blue 
dashed line: prediction for the critical radius at which cracks appear tp jR— ^/^eq|E q f (Eq j^eq) 
(equation (4.2). The solid (resp. empty) symbols correspond to the experimental data from 
hgure 2 (resp. figure 3), from which we determine the elastic moduli Eo{uo) (resp. (cracks)). 
Same symbols and colors as in figure 2a. 


4.1. Deviation to Taylor-Culick law 

The dynamics of the rim is then controlled by the balance between inertia and surface 
tension spatial gradient. We assume here that surface elasticity is constant up to a certain 
compression. In this particular case, the velocity of the aureole front (delimiting the 
frontier with the zone of undisturbed film whose thickness is still h = ho) simply reads 
Uf = y/2Eo/(pho) = Vc\/^o/7eq: which can be seen as a two dimensional analogous of 
sound (compression) velocity. The opening hole velocity can also be determined by solving 
the self-similar profile of the aureole and applying mass conservation. No analytical 
solution is provided in the considered radial geometry but numerical resolution shows 
that 

Uo = Vc f{Eo/jeq) ( 4 . 1 ) 

with / a decreasing function determined numerically (see appendix A) and reported in 
figure 4. It is thus still proportional to Taylor-Culick velocity Vc and decreases with the 
interfacial elasticity Eq, which is consistent with experimental observations of figure 2. 
From these data and equation (4.1), an interfacial elasticity Eo{uo) can be deduced 
(figure 4), which is reported in table 1 as a function of the MAc concentration. These 
data are compared to measurements of surface moduli F^od from the oscillating drop 
method performed by Mitrinova et al. (2013a). It shows the same qualitative variation 
with C despite a discrepancy on the absolute values obtained. However, the shrinkage 
amplitude and the compression timescales differ by several orders of magnitude, and the 
surfactant monolayer at the interface is expected to be highly non-Newtonian (Costa 
2012; Lucassen & Van Den Tempel 1972). 

4.2. Crack appearance 

Besides, snapshot inspection shows that cracks appear when the compressive surface 
wave (i.e. the aureole front) reaches the metallic frame of the film. Cracks are thus 
expected for: 



(4.2) 
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This prediction, represented in figure 4, is indeed in good agreement with our obser¬ 
vations: The hole radius when cracks appear Vp increases with the frame radius R and 
decreases with surface elasticity probed through MAc concentration variations, as shown 
in figure 3. Eventually, this critical compression does not depend on film thickness ho, 
showing that elasticity is not affected by confinement in the experimental configuration 
tested. 

The surface elasticity Eq (cracks) can therefore be deduced from the critical radius 
for crack apparition Vp (figure 4), and reported for the different MAc concentration in 
table 1. In this case, a very good agreement is obtained with the measured value of the 
surface modulus (Mitrinova et al 2013a), which confirms that the cracks arise from a 
compression of the aureole when its front reaches the frame. 

Note that the values of surface elasticity deduced from our two methods may differ. 
This is however expected due to our strong hypothesis of constant elasticity. Indeed, while 
the aureole front velocity only depends on the surface elasticity at very low compression 
rate (at the edge of the undisturbed film), the hole opening velocity modeling takes into 
account the elasticity through large interface compression. For large deformation, it is 
expected that the constant elasticity model fails: at large compression, the myristic acid 
surface concentration increases, which should result in larger elasticity as can be inferred 
from the moduli dependency upon C (Mitrinova et al. 2013a). The effective modulus 
Eq{uq) should then deviate more from measurements at small deformations performed 
by the oscillating bubble technique (Mitrinova et al. 2013a). 

In addition, the effect of elasticity has indirect consequences of some features of foam 
film rupture. For example, no flapping nor transverse destabilization of the rim was 
observed for our rigid soap films, in contrast to observations on low elasticity films 
and theoretical predictions (Lhuissier & Villermaux 20096); however, the reduced rim 
velocity could prevent the flapping instability to develop and subsequent film atomization 
(Lhuissier & Villermaux 20096). 


4.3. Cmek-like patterns 

Let us now discuss the observed crack-like patterns. During the fast deformation of 
the surface, the surfactants behave as an insoluble monolayer, comparable to a lipid 
monolayer experiencing a compression in a Langmuir trough (Lee 2008). In this case, 
above a critical compression, such a monolayer can behave differently depending on its 
structure. If it is liquid-like, it ejects the molecules in the bulk in the form of vesicles or 
bilayers. If it is solid-like, it can bend as an elastic sheet or fracture as a fragile material. 

Although our experiment does not provide a microscopic characterization of this 
transient surface structure, the crack pattern can be macroscopically characterized. In 
particular, even though the cracks are irregularly distributed, a number of cracks per 
radial segment can be counted; the deduced characteristic length between two cracks 
denoted A (figure Ic) is reported in figure 5 as a function of MAc concentration C (a) 
and film thickness ho (b). 

The increase of A with C is expected whatever the mechanism proposed. On the 
one hand, for higher bulk concentration, solubilization of interfacial surfactants is more 
difficult, hence a reduced number of vesicles or bilayers are to be expelled. On the other 
hand, a more concentrated solid-like layer will also exhibit a higher bending modulus 
and wavelength of elastic ripples are expected to increase with this modulus (Cerda & 
Mahadevan 2003). The decrease of the characteristic length with the film thickness ho 
is more unexpected. For the solid-like behavior, a thinner elastic sheet will bend more 
easily than a thicker one, thus exhibiting smaller ripple wavelength when buckled (Landau 
& Lifshitz 1975), in contrast with our observations. If cracks correspond to monolayer 
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Figure 5. (a) Characteristic length A between two cracks as a function of MAc concentration 
C for ho = 11 ±3 /am. (b) A as a function of the film thickness ho for solution D. (c) Bursting of 
a soap film of thickness ho = 3 /am (solution E). The timelapse between the two images is 9 ms, 
and the two lines highlight the velocity inhomogeneities. 


collapse by vesicles formation, it should not be affected by the film thickness. However, 
when modifying the thickness of the film, we also vary the velocity of compression or 
shrinkage rate. This parameter induces dynamical structural change in the surfactant 
monolayers (as it does in bulk crystallization processes for example (Cabane & Henon 
2003)). 

Finally, a complete understanding of the origin of these crack-like patterns would 
require some local high speed imaging structural analysis, which are beyond the scope 
of the present paper. 

The presence of these irregular cracks have direct consequences on hole opening 
dynamics. Indeed, when the aureole reaches the metallic frame, the hole opening slows 
down (and even stops for the thinner rigid films) and then irregularly accelerates in the 
region where the cracks appears. This feature is reported in figure 5c. Moreover, a velocity 
discontinuity in the liquid is observed, the outer region being at rest whereas the inner 
region is deformed. 


5. Conclusion 

To conclude, we have shown that modifying the chemistry of surfactant solutions 
can have strong influences on macroscopic dynamical processes, as observed in various 
situations in foams and foam films (Couder et al. 1989; Durand & Stone 2006; Seiwert 
ct al. 2013; Petit et al. 2015; Lorenceau et al. 2009). However, we have investigated 
here this effect under large deformations and in a fast dynamical process, i.e. at large 
Reynolds numbers, where the effects of molecular scales and surfactants are expected to 
be negligible. 

The initial constant velocity opening dynamics is well described taking into account the 
surface elasticity of the interfaces and was shown to be reduced at high surface modulus. 
This may be responsible for the inhibition of rim fragmentation and droplet ejection 
usually reported in liquid film ruptures (Lhuissier & Villermaux 20096). Further studies 
should determine the role of the ejected droplets in rupture propagation in macroscopic 
foams; the stability of these systems is indeed known to depend dramatically on the 
surface elastic properties (Rio & Biance 2014). However, finite size effects becomes soon 
crucial: when the elastic compression surface wave reaches the border of the frame, 
crack-like patterns, where velocity discontinuity are observed, appear in the foam film. 
Determining the origin of cracks, their microscopic structure, their location and number, 
and how they control film opening dynamics remain a challenge to tackle. 

The authors thank Gilles Simon for his help in setting up the experiment. 





P. C. Petit et al 



Figure 6 . (a) Profile of the film and notations, (b) Variations of surface tension 7 
versus shrinkage a in the simplihed constant elasticity modeling, (c) h/ho as a function of 
T/iVct) = y/2wlVc foi* radial bursting and different values of ac (0.2, 0.4, 0.6, 0.8 from top to 
bottom at the origin). 


Appendix A 

Equations for radial bursting 

We describe the radial bursting dynamics of a foam film of initial uniform thickness 
ho and include the effect of dynamic surface tension as first proposed by Frankel & 
Mysels (1969): the surface tension 7 is assumed to depend only on the shrinkage of the 
surface a which by mass conservation is related to film thickness a = ho/h. We denote 
the surface elasticity E{a) = ^. As viscous terms are negligible, the capillary forces 
are balanced by the fluid inertia. Variations of fluid velocity across the film are also 
neglected and equations are averaged over h. These equations can be explicitly solved in 
the unidimensional case (Frankel & Mysels 1969). However, in the case of radial bursting, 
a numerical resolution is necessary. 

We consider a material element that has initially the position R {i.e. that has La- 
grangian variables {RR)) . At instant t, its position is r{RR) and its thickness h{RR) 
(figure 6 a). The fluid velocity is u = drjdt and the shrinkage is defined as g = 
h/hQdr‘^/dR^ = {r/R)dr/dR. The momentum balance on the fluid element yields 


du d-f 

a = 


which can be rewritten 

du 2E{a) r da 2 ^ 

m ^ pho rM ^^^rM 

in which we have defined the characteristic velocity 


(Al) 


Ua = 



Following the analysis of Frankel & Mysels (1969), we are looking for self-similar 
solutions in the form r/t = f{R/t). We define the variables W = R‘^l{2t‘^) and w = 
r^/( 2 t^) {w and W have the dimensions of square velocities) and we expect w = w{W). 
The relative shrinkage is also set hy a = dw/dW. Starting from equation (Al), we find 


w 

W dw~ 

dw 


d‘^w 

w 


dW ^ 

^a=dw/dW ^ 

dW^ 


A first information on the film dynamics can be inferred from this equation: Far 
from the hole, i.e. for large VF, the film should remain undisturbed, which corresponds 
to w = W and dw/dW = 1. This condition combined with equation (A 2 ) yields 
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^a=i “ ~ which implies that the matching with the disturbed film can 

only be done at Vb = Wq = U‘f^^i/2. The velocity of the front of the aureole^ or extended 
rim corresponding to the disturbed film, is thus given hy Uf = Ua=i (Frankel & Mysels 
1969). 

Finally, the complete aureole profile and hole receding velocity will depend on the form 
of the elasticity versus shrinkage. 


Numerical resolution for a constant elasticity model 

We consider at first order a model of constant elasticity as described in figure 6b. 
We introduce here ac, which corresponds to the maximum shrinkage the film can endorse. 
For a > Q^c, the surface elasticity is constant and reads 


dy 

da 


— -Fq — Teq 


1 

1-ac 


and then 



where W = i/27eq/ (pho) is Culick velocity. 

When a > ac, equation (A 2) can be written as 


W r W drc 
w w dW 


dw 

dW 



2W^' 

w 


d^w 

dll^ 


(A3) 


In non-dimensionalized form (stating W = 2W/V^ and w = 2wjV‘f;)^ this equation 
reduces to 


w 

W dw' 

dtZi 


~ 


w 




2W^ d^w 


(A 4) 


with the two following boundary conditions. In VF = 0, at the hole, we have the maximum 
shrinkage (minimum value of ac): = 0) = ac. In W = Wq = 2 (i-a ) ^ 

aureole front, the solution should match the undisturbed film solution w{Wo) = Wq. 

This equation is solved numerically with a shooting method. From the function w{W), 
we can deduce the thickness profile, using the relation h/ho = l/{dw/dW) for different 
elasticities (figure 6c). As the elasticity increases (i.e. as ac becomes closer to I), we find 
that the aureole is thinner and wider, while for Eq = 7eq (corresponding to o^c = 0), one 
recovers a punctual rim receding at Taylor-Culick velocity 14- 

We can also estimate the initial hole velocity uq = ^/2w{W = 0), which is shown in 
figure 4 as a function of the ratio F^o/7eq- We also observe that the results obtained deviate 
from those obtained for unidimensional bursting (Frankel & Mysels 1969), especially for 
large elasticities, emphasizing the crucial role of radial geometry. 


REFERENCES 

Bird, James C., de Ruiter, Rielle, Courbin, Laurent & Stone, Howard A. 2010 
Daughter bubble cascades produced by folding of ruptured thin films. Nature 465 (7299), 
759-762. 

Cabane, B. & Henon, S. 2003 Liquides- Solutions, dispersions, gels. Belin. 

Cerda, E. & Mahadevan, L. 2003 Geometry and physics of wrinkling. Physieal Review Letters 
90 (7), 074302. 



10 P. C. Petit et al 

Costa, Severine 2012 Rheologie multi-echelle des mousse liquides. PhD thesis, Universite 
Paris-Est Marne La Vallee. 

COUDER, Y., Chomaz, J. M. & Rabaud, M. 1989 On the hydrodynamics of soap films. Phys. 
D 37 (1-3), 384-405. 

CULICK, F. E. C. 1960 Comments on a ruptured soap film. Journal of Applied Physics 31 (6), 
1128-1129. 

Durand, M. & Stone, H. A. 2006 Relaxation time of the topological T1 process in a two- 
dimensional foam. Physical Review Letters 97, 226101. 

Florence, A. T. & Frens, G. 1972 Aureole profile in bursting soap films - Surface-tension 
and surface relaxation in rapidly compressed monolayers. Journal of Physical Chemistry 
76 (21), 3024-3029. 

Florence, Alexander T. & Mysels, Karol J. 1974 Bursting of soap films. VI. effect of 
surfactant purity. The Journal of Physical Chemistry 78 (3), 234-235. 

Frankel, S. & Mysels, K. J. 1969 Bursting of soap films .2. Theoretical considerations. 
Journal of Physical Chemistry 73 (9), 3028-3038. 

CoLEMANOV, K., Denkov, N. D., Tcholakova, S., Vethamuthu, M. & Lips, A. 2008 
Surfactant mixtures for control of bubble surface mobility in foam studies. Langmuir 
24 (18), 9956-9961. 

Landau, L. D. & Lieshitz, E. M. 1975 Elasticity Theory. Oxford: Pergamon Press. 

Lastakowski, H., Boyer, F., Biance, A. L., Pirat, C. & Ybert, C. 2014 Bridging local to 
global dynamics of drop impact onto solid substrates. Journal of Fluid Mechanics 747, 
103-118. 

Lee, K. Y. C. 2008 Collapse mechanisms of Langmuir monolayers. Annu. Rev. Phys. Chem. 
59, 771-791. 

Lhuissier, H. & ViLLERMAUX, E. 2009a Destabilization of flapping sheets: The surprising 
analogue of soap films. Comptes Rendus Mecanique 337 (6-7), 469-480. 

Lhuissier, H. & Villermaux, E. 20096 Soap films burst like flapping flags. Physical Review 
Letters 103 (5), 054501. 

Liang, N. Y., Chan, C. K. & Choi, H. J. 1996 Dynamics of the formation of an aureole in 
the bursting of soap films. Physical Review E 54 (4), R3117-R3120. 

Lorenceau, E., Louvet, N., Rouyer, F. & Pitois, O. 2009 Permeability of aqueous foams. 
Eur. Phys. J. E 28, 293-304. 

Lucassen, j. & Van Den Tempel, M. 1972 Dynamic measurements of dilational properties 
of a liquid interface. Chemical Engineering Science 27 (6), 1283 - 1291. 

McEntee, W. R. & Mysels, K. J. 1969 Bursting of soap films .i. An experimental study. 
Journal of Physical Chemistry 73 (9), 3018-3028. 

Mitrinova, Z., Tcholakova, S., Golemanov, K., Denkov, N., Vethamuthu, M. & 
Ananthapadmanabhan, K. P. 2013a Surface and foam properties of SEES plus CAPB 
plus fatty acid mixtures: Effect of pH for C12-C16 acids. Colloids and Surfaces A- 
physicochemical and Engineering Aspects 438, 186-198. 

Mitrinova, Z., Tcholakova, S., Popova, Z., Denkov, N., Dasgupta, Bivash R. & 
Ananthapadmanabhan, K. P. 20136 Efficient control of the rheological and surface 
properties of surfactant solutions containing C8-C18 fatty acids as cosurfactants. Langmuir 
29 (26), 8255-8265. 

Petit, P., Seiwert, J., Cantat, I. & Biance, A.-L. 2015 On the generation of a foam film 
during a topological rearrangement. Journal of Fluid Mechanics 763, 286-301. 

Rio, Emmanuelle V Biance, Anne-Laure 2014 Thermodynamic and mechanical timescales 
involved in foam film rupture and liquid foam coalescence. ChemPhysChem 15 (17), 3692- 
3707. 

Seiwert, J., Monloubou, M., Dollet, B. & Cantat, I. 2013 Extension of a suspended soap 
film: A homogeneous dilatation followed by new film extraction. Physical Review Letters 
111 (9), 094501. 

Taylor, G. 1959 The dynamics of thin sheets of fluid .3. Disintegration of fluid sheets. 
Proceedings of the Royal Society of London Series A-Mathematical and Physical Sciences 
253 (1274), 313-321. 



